Timeline of changes in spike conformational dynamics in emergent SARS-CoV-2 variants reveal progressive stabilization of trimer stalk with altered NTD dynamics.

SARS-CoV-2 emergent variants are characterized by increased viral fitness and each shows multiple mutations predominantly localized to the spike (S) protein. Here, amide hydrogen/deuterium exchange mass spectrometry has been applied to track changes in S dynamics from multiple SARS-CoV-2 variants. Our results highlight large differences across variants at two loci with impacts on S dynamics and stability. A significant enhancement in stabilization first occurred with the emergence of D614G S followed by smaller, progressive stabilization in subsequent variants. Stabilization preceded altered dynamics in the N-terminal domain, wherein Omicron BA.1 S showed the largest magnitude increases relative to other preceding variants. Changes in stabilization and dynamics resulting from S mutations detail the evolutionary trajectory of S in emerging variants. These carry major implications for SARS-CoV-2 viral fitness and offer new insights into variant-specific therapeutic development.

Pharmacotherapeutic Potential of Natural Products to Target the SARS-CoV-2 PLpro Using Molecular Screening and Simulation Approaches.

Because of the essential role of PLpro in the regulation of replication and dysregulation of the host immune sensing, it is considered a therapeutic target for novel drug development. To reduce the risk of immune evasion and vaccine effectiveness, small molecular therapeutics are the best complementary approach. Hence, we used a structure-based drug-designing approach to identify potential small molecular inhibitors for PLpro of SARS-CoV-2. Initial scoring and re-scoring of the best hits revealed that three compounds NPC320891 (2,2-Dihydroxyindene-1,3-Dione), NPC474594 (Isonarciclasine), and NPC474595 (7-Deoxyisonarciclasine) exhibit higher docking scores than the control GRL0617. Investigation of the binding modes revealed that alongside the essential contacts, i.e., Asp164, Glu167, Tyr264, and Gln269, these molecules also target Lys157 and Tyr268 residues in the active site. Moreover, molecular simulation demonstrated that the reported top hits also possess stable dynamics and structural packing. Furthermore, the residues' flexibility revealed that all the complexes demonstrated higher flexibility in the regions 120-140, 160-180, and 205-215. The 120-140 and 160-180 lie in the finger region of PLpro, which may open/close during the simulation to cover the active site and push the ligand inside. In addition, the total binding free energy was reported to be - 32.65 ± 0.17 kcal/mol for the GRL0617-PLpro, for the NPC320891-PLpro complex, the TBE was - 35.58 ± 0.14 kcal/mol, for the NPC474594-PLpro, the TBE was - 43.72 ± 0.22 kcal/mol, while for NPC474595-PLpro complex, the TBE was calculated to be - 41.61 ± 0.20 kcal/mol, respectively. Clustering of the protein's motion and FEL further revealed that in NPC474594 and NPC474595 complexes, the drug was seen to have moved inside the binding cavity along with the loop in the palm region harboring the catalytic triad, thus justifying the higher binding of these two molecules particularly. In conclusion, the overall results reflect favorable binding of the identified hits strongly than the control drug, thus demanding in vitro and in vivo validation for clinical purposes.

SARS-CoV-2 Variants Increase Kinetic Stability of Open Spike Conformations as an Evolutionary Strategy.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) harbor mutations in the spike (S) glycoprotein that confer more efficient transmission and dampen the efficacy of COVID-19 vaccines and antibody therapies. S mediates virus entry and is the primary target for antibody responses, with structural studies of soluble S variants revealing an increased propensity toward conformations accessible to the human angiotensin-converting enzyme 2 (hACE2) receptor. However, real-time observations of conformational dynamics that govern the structural equilibriums of the S variants have been lacking. Here, we report single-molecule Förster resonance energy transfer (smFRET) studies of critical mutations observed in VOCs, including D614G and E484K, in the context of virus particles. Investigated variants predominately occupied more open hACE2-accessible conformations, agreeing with previous structures of soluble trimers. Additionally, these S variants exhibited slower transitions in hACE2-accessible/bound states. Our finding of increased S kinetic stability in the open conformation provides a new perspective on SARS-CoV-2 adaptation to the human population. IMPORTANCE SARS-CoV-2 surface S glycoprotein-the target of antibodies and vaccines-is responsible for binding to the cellular receptor hACE2. The interactions between S and hACE2 trigger structural rearrangements of S from closed to open conformations prerequisite for virus entry. Under the selection pressure imposed by adaptation to the human host and increasing vaccinations and convalescent patients, SARS-CoV-2 is evolving and has adopted numerous mutations on S variants. These promote virus spreading and immune evasion, partially by increasing the propensity of S to adopt receptor-binding competent open conformations. Here, we determined a time dimension, using smFRET to delineate the temporal prevalence of distinct structures of S in the context of virus particles. We present the first experimental evidence of decelerated transition dynamics from the open state, revealing increased stability of S open conformations to be part of the SARS-CoV-2 adaption strategies.

Integrating Conformational Dynamics and Perturbation-Based Network Modeling for Mutational Profiling of Binding and Allostery in the SARS-CoV-2 Spike Variant Complexes with Antibodies: Balancing Local and Global Determinants of Mutational Escape Mechanisms.

In this study, we combined all-atom MD simulations, the ensemble-based mutational scanning of protein stability and binding, and perturbation-based network profiling of allosteric interactions in the SARS-CoV-2 spike complexes with a panel of cross-reactive and ultra-potent single antibodies (B1-182.1 and A23-58.1) as well as antibody combinations (A19-61.1/B1-182.1 and A19-46.1/B1-182.1). Using this approach, we quantify the local and global effects of mutations in the complexes, identify protein stability centers, characterize binding energy hotspots, and predict the allosteric control points of long-range interactions and communications. Conformational dynamics and distance fluctuation analysis revealed the antibody-specific signatures of protein stability and flexibility of the spike complexes that can affect the pattern of mutational escape. A network-based perturbation approach for mutational profiling of allosteric residue potentials revealed how antibody binding can modulate allosteric interactions and identified allosteric control points that can form vulnerable sites for mutational escape. The results show that the protein stability and binding energetics of the SARS-CoV-2 spike complexes with the panel of ultrapotent antibodies are tolerant to the effect of Omicron mutations, which may be related to their neutralization efficiency. By employing an integrated analysis of conformational dynamics, binding energetics, and allosteric interactions, we found that the antibodies that neutralize the Omicron spike variant mediate the dominant binding energy hotpots in the conserved stability centers and allosteric control points in which mutations may be restricted by the requirements of the protein folding stability and binding to the host receptor. This study suggested a mechanism in which the patterns of escape mutants for the ultrapotent antibodies may not be solely determined by the binding interaction changes but are associated with the balance and tradeoffs of multiple local and global factors, including protein stability, binding affinity, and long-range interactions.

Dynamic Ca2+ sensitivity stimulates the evolved SARS-CoV-2 spike strain-mediated membrane fusion for enhanced entry.

Mutations in the spike protein generated a highly infectious and transmissible D614G variant, which is present in newly evolved fast-spreading variants. The D614G, Alpha, Beta, and Delta spike variants of SARS-CoV-2 appear to expedite membrane fusion process for entry, but the mechanism of spike-mediated fusion is unknown. Here, we reconstituted an in vitro pseudovirus-liposome fusion reaction and report that SARS-CoV-2 wild-type spike is a dynamic Ca2+ sensor, and D614G mutation enhances dynamic calcium sensitivity of spike protein for facilitating membrane fusion. This dynamic calcium sensitivity for fusion is found to be higher in Alpha and Beta variants and highest in Delta spike variant. We find that efficient fusion is dependent on Ca2+ concentration at low pH, and the fusion activity of spike dropped as the Ca2+ level rose beyond physiological levels. Thus, evolved spike variants may control the high fusion probability for entry by increasing Ca2+ sensing ability.

Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability.

Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex.

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

Microsecond simulations and CD spectroscopy reveals the intrinsically disordered nature of SARS-CoV-2 spike-C-terminal cytoplasmic tail (residues 1242-1273) in isolation.

All available SARS-CoV-2 spike protein crystal and cryo-EM structures have shown missing electron densities for cytosolic C-terminal regions (CTR). Generally, the missing electron densities point towards the intrinsically disordered nature of the protein region (IDPR). This curiosity has led us to investigate the cytosolic CTR of the spike glycoprotein of SARS-CoV-2 in isolation. The spike CTR is supposed to be from 1235 to 1273 residues or 1242-1273 residues based on our used prediction. Therefore, we have demonstrated the structural conformation of cytosolic region and its dynamics through computer simulations up to microsecond timescale using OPLS and CHARMM forcefields. The simulations have revealed the unstructured conformation of cytosolic region. Further, we have validated our computational observations with circular dichroism (CD) spectroscopy-based experiments and found its signature spectra at 198 nm. We believe that our findings will surely help in understanding the structure-function relationship of the spike protein's cytosolic region.

Investigating the conformational dynamics of SARS-CoV-2 NSP6 protein with emphasis on non-transmembrane 91-112 & 231-290 regions.

The NSP6 protein of SARS-CoV-2 is a transmembrane protein, with some regions lying outside the membrane. Besides a brief role of NSP6 in autophagosome formation, this is not studied significantly. Also, there is no structural information available to date. Based on the prediction by TMHMM server for transmembrane prediction, it is found that the N-terminal residues (1-11), middle region residues (91-112), and C-terminal residues (231-290) lies outside the membrane. Molecular Dynamics (MD) simulations showed that NSP6 consists of helical structures. In contrast, the membrane outside lying region (91-112) showed partial helicity, which was further used as a model and obtained disordered type conformation during 1.5 µs. Additionally, a 200ns simulation study of residues 231-290 have shown significant conformational changes. As compared to helical and beta-sheet conformations in its structure model, the 200ns simulation resulted in the loss of beta-sheet structures while helical regions remained intact. Further, we have experimentally characterized the residue 91-112 by using reductionist approaches. CD spectroscopy suggests that the NSP6 (91-112) is disordered-like region in isolation, which gains helical conformation in different biological mimic environmental conditions. These studies can be helpful to study NSP6 (91-112) interactions with host proteins, where different protein conformations might play a significant role. The present study adds up more information about the NSP6 protein aspect, which could be exploited for its host protein interaction and pathogenesis.

SARS-CoV-2 NSP1 C-terminal (residues 131-180) is an intrinsically disordered region in isolation.

The NSP1- C terminal structure in complex with ribosome using cryo-EM is available now, and the N-terminal region structure in isolation is also deciphered in literature. However, as a reductionist approach, the conformation of NSP1- C terminal region (NSP1-CTR; amino acids 131-180) has not been studied in isolation. We found that NSP1-CTR conformation is disordered in an aqueous solution. Further, we examined the conformational propensity towards alpha-helical structure using trifluoroethanol, we observed induction of helical structure conformation using CD spectroscopy. Additionally, in SDS, NSP1-CTR shows a conformational change from disordered to ordered, possibly gaining alpha-helix in part. But in the presence of neutral lipid DOPC, a slight change in conformation is observed, which implies the possible role of hydrophobic interaction and electrostatic interaction on the conformational changes of NSP1. Fluorescence-based studies have shown a blue shift and fluorescence quenching in the presence of SDS, TFE, and lipid vesicles. In agreement with these results, fluorescence lifetime and fluorescence anisotropy decay suggest a change in conformational dynamics. The zeta potential studies further validated that the conformational dynamics are primarily because of hydrophobic interaction. These experimental studies were complemented through Molecular Dynamics (MD) simulations, which have shown a good correlation and testifies our experiments. We believe that the intrinsically disordered nature of the NSP1-CTR will have implications for enhanced molecular recognition feature properties of this IDR, which may add disorder to order transition and disorder-based binding promiscuity with its interacting proteins.

Single-Molecule FRET Imaging of Virus Spike-Host Interactions.

As a major surface glycoprotein of enveloped viruses, the virus spike protein is a primary target for vaccines and anti-viral treatments. Current vaccines aiming at controlling the COVID-19 pandemic are mostly directed against the SARS-CoV-2 spike protein. To promote virus entry and facilitate immune evasion, spikes must be dynamic. Interactions with host receptors and coreceptors trigger a cascade of conformational changes/structural rearrangements in spikes, which bring virus and host membranes in proximity for membrane fusion required for virus entry. Spike-mediated viral membrane fusion is a dynamic, multi-step process, and understanding the structure-function-dynamics paradigm of virus spikes is essential to elucidate viral membrane fusion, with the ultimate goal of interventions. However, our understanding of this process primarily relies on individual structural snapshots of endpoints. How these endpoints are connected in a time-resolved manner, and the order and frequency of conformational events underlying virus entry, remain largely elusive. Single-molecule Förster resonance energy transfer (smFRET) has provided a powerful platform to connect structure-function in motion, revealing dynamic aspects of spikes for several viruses: SARS-CoV-2, HIV-1, influenza, and Ebola. This review focuses on how smFRET imaging has advanced our understanding of virus spikes' dynamic nature, receptor-binding events, and mechanism of antibody neutralization, thereby informing therapeutic interventions.

Surveying the Side-Chain Network Approach to Protein Structure and Dynamics: The SARS-CoV-2 Spike Protein as an Illustrative Case.

Network theory-based approaches provide valuable insights into the variations in global structural connectivity between different dynamical states of proteins. Our objective is to review network-based analyses to elucidate such variations, especially in the context of subtle conformational changes. We present technical details of the construction and analyses of protein structure networks, encompassing both the non-covalent connectivity and dynamics. We examine the selection of optimal criteria for connectivity based on the physical concept of percolation. We highlight the advantages of using side-chain-based network metrics in contrast to backbone measurements. As an illustrative example, we apply the described network approach to investigate the global conformational changes between the closed and partially open states of the SARS-CoV-2 spike protein. These conformational changes in the spike protein is crucial for coronavirus entry and fusion into human cells. Our analysis reveals global structural reorientations between the two states of the spike protein despite small changes between the two states at the backbone level. We also observe some differences at strategic locations in the structures, correlating with their functions, asserting the advantages of the side-chain network analysis. Finally, we present a view of allostery as a subtle synergistic-global change between the ligand and the receptor, the incorporation of which would enhance drug design strategies.

Real-Time Conformational Dynamics of SARS-CoV-2 Spikes on Virus Particles.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) mediates viral entry into cells and is critical for vaccine development against coronavirus disease 2019 (COVID-19). Structural studies have revealed distinct conformations of S, but real-time information that connects these structures is lacking. Here we apply single-molecule fluorescence (Förster) resonance energy transfer (smFRET) imaging to observe conformational dynamics of S on virus particles. Virus-associated S dynamically samples at least four distinct conformational states. In response to human receptor angiotensin-converting enzyme 2 (hACE2), S opens sequentially into the hACE2-bound S conformation through at least one on-path intermediate. Conformational preferences observed upon exposure to convalescent plasma or antibodies suggest mechanisms of neutralization involving either competition with hACE2 for binding to the receptor-binding domain (RBD) or allosteric interference with conformational changes required for entry. Our findings inform on mechanisms of S recognition and conformations for immunogen design.